Spontaneous Unknotting of a Polymer Confined in a Nanochannel 
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We study the dynamics of a knot in a semiflexible polymer confined to a narrow channel of 
width comparable to the polymers' persistence length. Using a combination of Brownian dynamics 
simulations and a coarse-grained stochastic model, we characterize the coupled dynamics of knot size 
variation and knot diffusion along the polymer, which ultimately leads to spontaneous unknotting. 
We find that the knot grows to macroscopic size before disappearing. Interestingly, an external force 
applied to the ends of the confined polymer speeds up spontaneous unknotting. 



With recent nanotechnology, single biopolymers can be 
manipulated and observed inside fabricated nanoscale de- 
vices Jj]. Such 'lab-on-a-chip' techniques are interesting 
to study the physical or biochemical properties of indi- 
vidual molecules, and to sort them according to these 
properties. A key benefit of confining, e.g., a long DNA 
molecule in a nanochannel is the resulting stretched-out 
conformation which permits the mapping of a position 
along the channel axis to a position on the polymer con- 
tour. Thereby, confinement offers the possibility to di- 
rectly observe and map local sequence-dependent prop- 
erties of DNA such as affinity for protein binding, which 
can conventionally be studied only indirectly with bio- 
chemical methods. The physics of polymers in confined 
geometries 0, 0, 0, Q is fundamental to these applica- 
tions. An interesting problem arises from topological de- 
fects in the form of knots in the polymer contour 0, 0] , 
which destroy the mapping between channel and poly- 
mer coordinates. Knots easily form in a linear polymer, 
before or while it is threaded into a nanochannel [8|. 



Here, we consider the situation after the polymer is 
trapped inside the channel. Since a knot in a linear 
polymer is not topologically conserved, it eventually dis- 
appears spontaneously, driven by thermal fluctuations 
0,11, E, U I ■ Our focus is on the physics of this process 
under strong confinement. The typical conformation of 
a linear polymer with contour length L inside a chan- 
nel strongly depends on the channel width d. It forms a 
random coil like a free polymer, if the channel is wider 
than the polymers' radius of gyration (and filled with a 
good solvent). In more narrow channels, the conforma- 
tion elongates as described within the blob picture [l^ . 
Further reducing d to the scale of the polymers' persis- 
tence length ip gradually stretches out the polymer, by 
constraining the local tangent vectors t(s) of its contour 
(0 < s < L) to directions close to the direction c of the 
channel axis [^. Concomitantly, the typical number of 
"U-turns" , i.e., sign reversals of t(s) • c, decreases to zero 
[3| • Recent experiments have entered this strong confine- 
ment regime with DNA molecules j^. In this regime, the 
characteristic lengthscale of the polymer conformation is 
the Odijk length Ld ^ (d'^lpY^^, the typical distance be- 



tween subsequent collisions of the polymer contour with 
the channel wall 0. 

In a knotted configuration, at least two U-turns are 
topologically interlocked, see Fig. [1] (top). These U-turns 
can disappear only at the polymer ends (their "pair anni- 
hilation" is impossible), either at the same or at different 
ends. In terms of the knot, these two options correspond 
to spontaneous unknotting via diffusion of the knot to 
one end, or to swelling of the knot to macroscopic sizes 
of order L. Knot diffusion was experimentally observed 
in DNA molecules stretched by an externally applied 
tensio n [71 , and addressed in recent theoretical studies 
[gl [loL [llf . However, the dynamics of confined knots was 
either not considered H, 11 1 or treated on an equal foot- 
ing with force- induced stretching We find that these 
two situations lead to qualitatively different dynamics. In 
particular, we identify a dynamic interplay between knot 
diffusion and swelling under strong confinement, which 
is absent for force-induced stretching. 

Due to the enormous range of timescales in the sys- 
tem, we take a multiscale approach to characterize its 
dynamics. The shortest relevant timescales are indepen- 
dent of L and correspond to the local relaxation of the 
polymer degrees of freedom in, e.g., a single Odijk length, 
whereas the longest relevant timescale is the lifetime of 
a knot, which increases strongly with the length of the 
polymer, see below. We use Brownian dynamics simu- 
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FIG. 1: (Top) Illustration of a knotted polymer confined to 
a nanochannel, and definition of the left, center, and right 
segment with lengths I, k, and r, respectively. (Bottom) Il- 
lustration of the three slow modes of motion allowing the knot 
to move and change its size. 
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lations to explore the full polymer dynamics over short 
periods, and then match the results to a coarse-grained 
stochastic model for the long term dynamics. 

To study the full polymer dynamics, we use a 'bead- 
spring' polymer model consisting of N beads at po- 
sitions Vi coupled by a harmonic spring potential. 
Us = KsY.ii\^i,i+i\ - ^)^/2 with = r,; - rj, 

which keeps neighboring beads approximately at a dis- 
tance b. A soft excluded volume potential, Uex — 
'ikBTJ2j<:i-i9{<^/\^i,j\) prevents self-crossing of the 
polymer. Here, a denotes the interaction range, and 
g{x) — (x^^ — 2x^ + 1)8(1 — x) yields the repulsive 
part of the Lennard- Jones potential with the Heav- 
iside function 9 (a;). The bending stiffness is de- 
scribed by the discretized worm-like chain energy Ub = 
ksTtp/b ■ ~ cos6'i) with the local bending angle 

9i at bead i. The polymer is confined to a channel 
with circular profile by a soft wall potential, Uch = 
ksTe-^ J2i (hi - ^+ e)*e{hi - | + e), where hi is the 
distance of bead i from the channel axis, d the channel 
diameter, and e parameterizes the softness of the poten- 
tial. We measure all lengths in units of d and energies 
in units of ksT, i.e., d = I and ksT = 1 in numerical 
simulations. To put us in the strong confinement regime, 
we choose a persistence length ip = 3 and consider only 
polymer lengths L = Nb ^ £p. Furthermore, we use a 
spring constant Kg — 10^, a segment length b = 0.1, a 
bead-bead interaction range a — 1.2b, and e — 0.05. 

Taken together, we have the energy function U = Us + 
Uex + Ub+Uch- The Brownian dynamics of the polymer is 
then described by the discrete-time Langevin equations, 



r,(t-f At) =r,(i)-/ibAiVr.L^+v/2Ai(rfcBTAtr7,(t) (1) 

with a bead mobility fib, a time step At, and random 
forces T]i{t) with a variance (r;j(t) • r]j{t')) — 36ijSu'- 
We choose our time unit such that a polymer segment of 
unit length has unit mobility, fib = b, and a time step 
At = 2 • 10~^, a tradeoff between accuracy and efficiency. 

Eq. Il]) describes the coupled dynamics of the beads, 
but we are ultimately interested in how the conformation 
evolves on a coarse-grained scale. To this end, we con- 
sider simple trefoil knots as in Fig. [Jl and monitor two 
observables, the knot size k and position p defined as 

fc EE 6^e(c • r,^j+i) , PEE — ^ie(c-rj,j+i) . (2) 



Here c is a vector along the channel axis, with the same 
orientation as r^va, see Fig. [1] These definitions are ap- 
plicable in the strong confinement regime, where spon- 
taneous formation of U-turns is suppressed and a trefoil 
knot has only two interlocked U-turns. The length k of 
the polymer contour between the turns is then a well- 
defined measure of the knot size, as is the location p of 
its midpoint for the knot position. 
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FIG. 2: Diffusion coefficients {Dp, Dk, Dkp, top) and drift 
velocities {vp,Vk, bottom) for various initial fco (left) and po 
(right) obtained from Brownian dynamics simulations (data 
points) 16^ and a coarse-grained model (lines). 



We use initial polymer configurations drawn from an 
ensemble with a specified knot position pq and size kg 
and evolve them according to Eq. We find 

that the knot diffuses along the polymer and at the same 
time changes its size, apparently in an erratic, non-biased 
way. To quantify this behavior, we consider the time- 
dependent position shift Sp(t) = p{t) — po and size shift 
6k(t) = k{t) — kQ, and determine their averages, variances, 
and covariance from 150 independent simulations for each 
combination of initial values fco and pq. For these simula- 
tions, we used total simulation times considerably longer 
than required for relaxation of the nearly straight poly- 
mer sections, but short enough to avoid large changes in 
p or fc within a simulation. The data is thus suitable 
to estimate local drift velocities Vp, Vk and diffusion co- 
efficients Dp, Dk in knot position and size space. The 
symbols in Fig. [2] show these estimates for a polymer of 
length L — 50, both as a function of fco at a fixed pq 
and as a function of po with fixed fco . Also shown is the 
cross-correlation coefficient D^p — dt{Sp6k) /2. 

The drift /diffusion coefficients Vp, Vk, Dp, Dk, and Dkp 
provide the "interface" between our detailed Brownian 
dynamics analysis and our coarse-grained description in 
(fc,p)-space. In the following we will see that the de- 
pendence on fco and po observed in Fig. [2] is indeed very 
plausible. We then study the implications of this depen- 
dence for the long term knot dynamics. 

On a larger scale, an equation for the stochastic dy- 
namics of the knot in (fc,p)-space can be derived if the 
relevant modes of motion and their associated mobilities 
are known. The long time dynamics is dominated by 
the slow modes where a large segment of the polymer 
moves relative to the solvent. For our trefoil knot, there 
are three segments, the central segment between the two 
U-turns with contour length fc, the left end with length 
I = p—k/2, and the right end with length r = L—p—k/2. 
As illustrated in Fig.[Tl each of these segments has an as- 
sociated mode: (i) In a "left mode", the left end moves 
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and exchanges length with the center segment, (ii) Simi- 
larly, in a "right mode" , the right end moves, exchanging 
length with the center, (iii) In contrast, a "knot mode" 
exchanges polymer length between the left and the right 
end, but only the center segment moves. The mobility of 
each mode scales inversely with the length of the segment 
moving relative to the solvent. 

The mode approach leads to coupled stochastic differ- 
ential equations for the knot position p(t) and size k{t) 
in the general Ito type form 



d 
di 



V D - M V U{k,p) + V2B r]{t) 



(3) 



Here, the diffusion matrix D is linked to the mobil- 
ity matrix M via the Einstein relation D — fcsTM. 
The potential U{k,p) can be used to study the effect 
of additional forces, see below, but is set to zero for 
now. The noise vector 77 has three uncorrelated com- 
ponents r]i, one for each mode, with (rjiit)) = and 
{'rii{t)rjj{t')) — Sij6{t — t'). The matrix B describes 
the effect of the three noise components on the coordi- 
nates {k,p). It is connected to the diffusion matrix via 
D = BB'^. Hence, knowledge of B suffices to specify the 
coarse-grained dynamics explicitly. For the above modes, 
we find 



B 



1 / 2/Vl 

4 V -i/Vi i/V^ 2/Vfc, 



(4) 



where we again set the mobility of a polymer stretch 
of unit length as well as ksT equal to one. The dif- 
fusion matrix and its derivatives yield simple analyti- 
cal expressions, without any adjustable parameters, for 
our drift/diffusion coefficients. Conversely, the five co- 
efficients fully specify the infinitesimal motion in {k,p)- 
space, and the long term dynamics is obtained by inte- 
grating Eq. ([3]) or the equivalent Fokker-Planck equation. 
Fig. [5] suggests that the three slow modes correctly cap- 
ture the coarse-grained dynamics of the knot (see, how- 
ever, the closer analysis below). 

The modes also facilitate the physical interpretation 
of Fig. [21 For instance, the knot mode dominates the 
mobility of the knot. Dp = Jq{^ + j + with a con- 
tribution that depends inversely on k (smaller knots are 
more mobile). Knot size fiuctuations increase when the 
knot approaches a polymer end, Dk = j{j + ^), and 
Dk also grows slowly with knot size due to the con- 
comitant shrinking of the ends. The cross-correlation 
Dkp = |(^ — j) is appreciable only close to the edges, 
where the dominant end mode simultaneously leads to a 
shift in the knot position and size. Within the mode ap- 
proach, drift coefficients arise only via the noise-induced 
drift term, VD. The size drift Vk — + :pi) tends 
to increase the knot size, whereas the positional drift 
Vp = — p-) tends to push knots to the edges. 



Since even a small drift velocity dominates over dif- 
fusion at long timescales, the knot dynamics is poten- 
tially very sensitive to any drift effects. Both drift effects 
discussed above are very small and our analysis below 
shows they are not relevant over the typical lifetime of 
a knot. However, in the simulation data in Fig. [21 the 
statistical errors for the drift coefficients are comparable 
to the average values, and both are much larger than the 
noise-induced drift. The positional drift values Wp(fco,Po) 
show no systematic bias and are compatible with zero or 
just the noise-induced drift. In contrast, the size drift 
values Vk{kQ,pQ) display a bias to negative values, sug- 
gesting systematic knot shrinkage. To pinpoint the origin 
of this effect, we first tested whether it depends on the 
topology of the polymer configuration: We repeated the 
simulations starting from initial configurations with two 
U-turns and the same A:o,po values, but without knot. 
These simulations resulted in the same bias as for the 
knotted configuations (data not shown). We then re- 
peated the simulations again, but with the excluded vol- 
ume potential Uex switched off. In this case, we observed 
no bias (data not shown). 

These results indicate that the self-exclusion between 
the three polymer strands in the channel region between 
the U-turns causes the bias. The self-exclusion effect is 
noticeable in our simulations due to the large bead size, 
which is required for computational efficiency. However, 
under typical experimental conditions it would be negli- 
gible except under extreme confinement. On the theoret- 
ical side, we expect no additional drift effects, since the 
entropy of the three polymer segments of Fig. [H is addi- 
tive f2] , and the bending energy of the U-turns does not 
depend on k. Note that the confinement is essential for 
these arguments; knots in unconfined semiflexible poly- 
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FIG. 3: (a) Distribution of final knot sizes K at unknotting, 
scaled by the polymer length L. The broad distribution shows 
that knots become macroscopically large, (b) Effect of an 
applied external force / on the mean knot lifetime (r) for 
L — 5000. Above a threshold the (r) becomes much smaller, 
i.e., unknotting is accelerated, (c) Scaling of mean final knot 
size (K) and mean knot lifetime (r) with polymer length L 
for / = and / = 0.05. 
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Next, we explore the long term knot dynamics. Since 
changes in k and p occur on similar timescales, see Fig. [2l 
we expect a dynamic interplay between these two degrees 
of freedom. Specifically, we address the question whether 
knots disappear by growing, by diffusing to the end, or 
by a combination of both processes. To this end, we em- 
ploy the coarse-grained eq. ([3]) and incorporate the phys- 
ical constraint that the knot cannot shrink to arbitrarily 
small size via the potential U(k,p) with a steeply rising 
barrier at fc = 1 . A natural definition of knot disappear- 
ance is min(^, r) < k. Starting with a small knot (fco = 5) 
in the center of the polymer {po = L/2), we determine the 
final knot size K just before unknotting. Fig. [3ja) shows 
the resulting distribution of K/L for different polymer 
lengths L. The broad range of final knot sizes indicates 
that polymers neither unknot by knot growth nor by knot 
diffusion alone, but by a combination of both processes, 
and that knots typically grow to macroscopic sizes ~ L 
before disappearing. 

This scaling behavior can be rationalized with a sim- 
ple argument. Consider a knot diffusing along the poly- 
mer until it dissolves after a time r. It reaches a typical 
size k ~ \/DkT with Dk on the order of 1/L. Hence, 
k"^ ^ t/L. On the other hand, the typical time until 
the knot reaches a polymer end scales as r ~ / Dp 
with Dp ~ 1/fc. Taken together, the knot has a char- 
acteristic final size k ^ L and lifetime t ^ . By the 
same arguments, one finds that the noise-induced drift 
is (marginally) irrelevant, i.e., it does not affect the scal- 
ing of r. Generally, any drift effect on the knot size is 
significant only when v], > Dk / r. 

How is the knot dynamics affected when a weak ex- 
ternal pulling force / is applied to the ends of the con- 
fined polymer? We consider the force regime where the 
knot size is constrained without tightening the knot so 
that molecular friction would become important. We 
include the pulling force through an additional term 
fir+l-k) = fi2k-L) in the potential C/(fc,p). Fig.Etb) 
shows the mean unknotting time (r) as a function of the 
applied force. While the unknotting time is independent 
of the applied force for small /, it then decreases, i.e., 
an increase in the tension accelerates unknotting. This 
result becomes intuitively clear when considering that 
forces much larger than ^ k^T/L constrain the knot to 
sizes much smaller than the polymer length; the knot 
then mainly diffuses along the polymer with constant 
mobility, leading to a typical unknotting time r X^, 
faster than the L'^-scaling without force, see Fig. [3ljc). 
For experimentally relevant parameters, the cross-over of 
Fig. [3jb) occurs at very small forces, e.g., in the sub-pN 
regime for a DNA of 1 /xm length. 

To summarize, we showed that knot dynamics under 
strong confinement is qualitatively different from knot 
dynamics under tension. While the polymer is stretched 
out in both cases, the knot remains localized in the lat- 
ter, whereas in the former, knot sizes and knot posi- 



tion change simultaneously, leading to knot growth up to 
macroscopic sizes. As a consequence, tiny forces which 
keep the knot localized can speed up unknotting consid- 
erably. Experimentally, such forces can easily occur, e.g., 
by weak attachment and hydrodynamic drag. 

Our theoretical treatment can be generalized to more 
complicated knots, which may be associated with more 
U-turns and associated modes. However, our model is 
clearly very simplistic in several respects. Most impor- 
tantly, we did not explicitly consider hydrodynamic and 
electrostatic interactions between the wall and the poly- 
mer, as well as between different parts of the polymer. 
These effects certainly influence the absolute timescale 
of the knot dynamics. However, we expect that they 
do not affect the dynamics qualitatively. For instance, 
electrostatic effects should be weak as long as the Debye 
screening leng th is considerably smaller than the channel 
diameter We hope that our assumptions and pre- 

dictions will be scrutinized experimentally and help to 
explore the physics of polymers in nanofluidic devices. 
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